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Probability density functions (PDFs) of scale-dependent energy fluctuations, P[5E(£)], are stud- 
ied in high-resolution direct numerical simulations of Navier-Stokes and incompressible magneto- 



hydrodynamic (MHD) turbulence. MHD flows with and without a strong mean magnetic field 

^J . are considered. For all three systems it is found that the PDFs of inertial range energy fluctu- 

ations exhibit self-similarity and monoscaling in agreement with recent solar-wind measurements 
O 
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[B. Hnat et al., Geophys. Res. Lett. 29(10), 86-1 (2002)]. Furthermore, the energy PDFs exhibit 
similarity over all scales of the turbulent system showing no substantial qualitative change of shape 
as the scale of the fluctuations varies. This is in contrast to the well-known behavior of PDFs of 
turbulent velocity fluctuations. In all three cases under consideration the P[5E(£)] resemble Levy- 
Z^ ■ type gamma distributions ~ A" 1 e~x.p(—\5E\/ A)\5E\~~' y The observed gamma distributions exhibit 

a scale-dependent width A(£) and a system-dependent 7. The monoscaling property reflects the 



inertial-range scaling of the Elsasser-field fluctuations due to lacking Galilei invariance of 5E. The 



appearance of Levy distributions is made plausible by a simple model of energy transfer. 
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Turbulence in electrically conducting magnetofluids is, apart from its importance for 
laboratory plasmas (see, for example 1]), a key ingredient in the dynamics of, e.g., the 
earth's liquid core and the solar wind (see e.g. 2|). A simple description of such plasmas 
is the framework of incompressible magnetohydrodynamics (MHD), a fluid approximation 
neglecting kinetic processes occuring on microscopic scales. This approach is appropriate if 
the main interest is focused on the nonlinear dynamics and the inherent statistical properties 
of fluid turbulence. To this end, two-point increments of a turbulent field component, say 
/, in the direction of a fixed unit vector e, Sf(£) = f(r + ei) — /(r) are analyzed, since they 
yield a comprehensive and scale-dependent characterization of the statistical properties of 
turbulent fluctuations via the associated probability density function (PDF) [3]. 

PDFs of temporal fluctuations [16] in the solar wind, e.g. of total (magnetic + kinetic) 
energy density, as measured by the WIND spacecraft are self-similar over all observed scales, 
exhibit monoscaling, and closely resemble gamma distributions. In contrast the PDFs of 
velocity and magnetic field are found to display well-known multifractal characteristics, i.e. 
the associated PDFs change from Gaussian at large scales to leptocurtic (fat-tailed) at small 
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scales |4|, |5|, 16]. The solar wind plasma is a complex and inhomogeneous mixture of mutu- 
ally interacting regions with different physical characteristics and dynamically important 
kinetic processes [7|, |8j. Thus, it is not clear if the abovementioned solar-wind observa- 
tions are caused by turbulence or some other physical phenomenon. This paper reports 
an investigation of turbulent PDFs based on high-resolution direct numerical simulations 
of physically 'simpler' homogeneous incompressible MHD and Navier-Stokes turbulence to 
elucidate this question. Monoscaling of the two-point PDFs of energy is found in the iner- 
tial range of macroscopically isotropic MHD turbulence, anisotropic MHD turbulence with 
an imposed mean magnetic field as well as in turbulent Navier-Stokes flow. The respective 
PDFs resemble leptocurtic gamma laws on all spatial scales in agreement with the solar- wind 
measurements. The monoscaling property is shown to be a consequence of lacking Galilei 
invariance of the energy fluctuations in combination with turbulent inertial-range scaling. 
The appearance of Levy-type gamma distributions apparently results from nonlinear turbu- 
lent tranfer as suggested by similar findings in all three investigated systems and a simple 
reaction-rate model. 

The dimensionless equations of incompressible MHD, formulated in Elsasser variables 
z^ = v ± b with the fluid velocity v and the magnetic field b which is given in Alfven-speed 



units 9], read 

V-z ± = (1) 

d t z ± = -z T ■ Vz ± - VP + j? + Az ± + ^Az T (2) 

with the total pressure P = p+ \b 2 . The dimensionless kinematic viscosity \x and magnetic 
diffusivity r\ appear in i]± = 1/2 (/i ± ??). 

The data used in this work stems from pseudospectral high-resolution direct numerical 
simulations [lCj based on a set of equations equivalent to Eqs. ([I]) and ([2]). It describes 
homogeneous fully-developed turbulent MHD and Navier-Stokes (b = 0) flows in a cubic 
box of linear size 2tc with periodic boundary conditions. The initial conditions for the 
decaying simulation run consist of random fluctuations with total energy equal to unity. 
In the MHD cases total kinetic and magnetic energy are approximately equal. The initial 
spectral energy distribution is peaked at small wavenumbers around k = 4 and decreases like 
a Gaussian towards small scales. In the MHD setups magnetic and cross helicity are small 
implying z + ~ z~ . The driven turbulence simulations were run towards quasi-stationary 
states whose energetic and helicity characteristics as mentioned above are roughly equal to 
the decaying run. The MHD magnetic Prandtl number Pr m = fi/r] is unity. The Reynolds 
numbers of all configurations are of order 10 3 . 

Three cases are considered. Setup (a) represents decaying macroscopically isotropic 3D 
MHD turbulence. The dataset contains 9 states of fully developed turbulence each compris- 
ing 1024 3 Fourier modes. The samples are taken equidistantly in time over a period of about 
3 large eddy turnover times. The angle-integrated energy spectrum of this system exhibits 



a Kolmogorov-like scaling law 11] in the inertial range, i.e. E/. ~ k~ 5 ^ 3 . The second dataset 
(b) contains simulation data of a driven quasi-stationary macroscopically anisotropic MHD 
flow with a strong constant mean magnetic field. The driving is accomplished by freezing 
the largest Fourier modes of the system (k < 2). The data comprises 1024 2 Fourier modes 
perpendicular to the direction of the mean field and 256 modes parallel to it. This dataset 
covers about 2 large eddy turnover times of quasi-stationary turbulence with 8 samples 
taken equidistantly over that period. The perpendicular energy spectrum shows Iroshnikov- 



Kraichnan-like behavior Eu. ~ k 
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131 ] . Note that this is neither claim nor clear 



evidence for the validity of the Iroshnikov-Kraichnan picture in this configuration. For fur- 

r 

ther details of the simulations and additional references see 1101. The third simulation (c) 



represents a turbulent statistically isotropic Navier-Stokes flow with resolution 1024 3 which 
is kept stationary by the same driving method as in case b) and exhibits Kolmogorov-scaling 
Ek ~ k~ 5 ^ 3 of the turbulent energy spectrum. 

For all turbulent systems the statistical properties of Sf(£) which is computed over vary- 
ing scale £ are investigated. In the present work / stands for the component of z + in 
the increment direction e or the fluctuation energy defined here as E = {z + ) 2 . For the 
macroscopically isotropic setups (a) and (c) e = e z . In case (b) the unit vector points in a 
fixed arbitrary direction perpendicular to the mean magnetic field. Under the assumption 
of statistical isotropy which is approximately fulfilled in setup (a), (c), and in system (b) in 
planes perpendicular to the mean magnetic field, the statistical properties of 5f(£) depend 
solely on £. This assumption also holds approximately for SE if contributions by eddies on 
larger-scales which are convolved into this non-Galileian-invariant quantity can be regarded 
as quasi- const ant on scale £ (see below). The quantity (8f(£)) scales self-similarly with the 
scaling parameter a (a > 0), if (Sf(X£)) = X a (f(£)) for every A. For the associated cumula- 
tive probability distribution follows p(Sf(£) < p) — p(X~ a 5f(X£) < p) for any real p. This 
implies for the probability density P 

P[5f(£)} = X- a P s [X- a 5f s ] (3) 

introducing the master PDF P s with 5f s = 5f(X£). According to Eq. ([3]), there is a family 
of PDFs that can be collapsed to a single curve P s , if a is independent of £. This is known 
as monoscaling in contrast to multifractal scaling observed, e.g., for two-point increments of 
a turbulent velocity field. 

To test if the abovementioned observations in the solar wind are a phenomenon related to 
inherent properties of turbulence time- and space-averaged increment series 5z + (£) and 5E(£) 
for different £, ranging between n/ 512 up to it are computed. In system (a) the increments 
are normalized using (E T y/ 2 with E T = 1/4 J v dV[(z + ) 2 + (z~) 2 } to compensate for the 
decaying amplitude of the turbulent fluctuations. The PDFs are generated as normalized 
histograms of the respective increments taken over all positions in the 27r-periodic box which 
contains the real space fields, v(r) and b(r), computed from the available Fourier-coefficients. 
Fig. [T]shows P[5E(£)] for various £ in the isotropic case (a). The non-Gaussian nature of the 
PDFs over all scales is evident. Similar behavior is found in the anisotropic case (b) where 
the increments are taken perpendicularly to the direction of the mean field as well as in 
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FIG. 1: The PDFs of total energy fluctuations 8E on five different scales £ = -n/n with n = 511 
(solid), n = 130 (dotted), n = 46 (dashed), n = 4 (dot-dashed), n = 1 (3-dot-dashed). 



the Navier-Stokes simulation (c). The PDFs are highly symmetric and become increasingly 
broader with growing £ reflecting the increase of turbulent energy towards largest scales. 
Interestingly, the PDFs at all scales have the same leptokurtic shape resembling Levy laws. 
In particular, away from the center, SE = 0, the PDFs are close to gamma distributions 
~ exp(— \SE\/A)\SE\~ J of different widths A. The exponent 7 of the best fits is constant in 
the inertial range and amounts approximately to 3.4 (a), 4.2 (b), and 3.1 (c). In the solar 
wind a similar finding however with 7 ^ 2.5 was reported 4J. 

The similarity of the P[SE{i)\ on different scales £ suggests the possibility of monoscaling. 
The monoscaling exponent is expected to be scale-independent in the inertial range only since 
the energy increments are not Galilei invariant. Therefore, small-scale SE also comprise 
contributions by larger eddies which advect the small-scale fluctuations. A linearization 
of SE with respect to the largest-scale contribution (z^) 2 ^> (Sz + ) 2 yields to lowest order 
SE « (zq + Sz + ) 2 ~ ZqSz + . As a consequence, the energy increments reflect the inertial- 
range scaling of the turbulent Elsasser fields, i.e. SE ~ Sz + ~ C a . To apply the rescaling 
procedure given by Eq. ([3]) (cf. also |4j) the exponent a is extracted from the PDFs by two 
independent techniques. 

Firstly, the standard deviation is considered which is defined as a(£) = [(SE(i) 2 )] 1 ^ 2 . In 
the inertial range a exhibits power-law behavior with respect to the increment distance, 
a(£) ~ £ a , Fig. [2] shows the standard deviation of total energy fluctuations in the inertial 
range for the isotropic case (a) in double logarithmic presentation. A linear least-squares 
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FIG. 2: Standard deviation of total energy increments within the inertial range in case (a) (trian- 
gles) with linear least-squares fit (solid line). 
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FIG. 3: Rescaled PDFs of total energy fluctuations in the inertial range of the isotropic case (a). 
The gamma law 10 -3 exp(— \5E\/0.35)\5E\^ S ' 1 is represented by the dashed curve. 

fit is carried out to obtain a. The characteristic exponents deduced in this way are a = 
0.29 ± 0.025 for the isotropic case (a), a = 0.23 ± 0.025 for the anisotropic case (b), and 
a = 0.28 ±0.03 for the Navier-Stokes flow (c). As expected these values are close to the non- 
intermittent scaling exponents observed for the turbulent field fluctuations, i.e. «K4i = 1/3 
for cases (a) and (c) while aiK = 1/4 for case (b). 

Secondly, in the inertial range the characteristic exponents can be obtained via the ampli- 
tude of P(0, £) ~ £~ a profiting from the fact that the peaks of the PDFs are statistically the 
least noisy part of the distributions. The scaling exponent obtained by using this method 
is in good agreement with the value of a obtained via the PDF variance. Fig. [3] shows 
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FIG. 4: Rescaled PDFs of total energy fluctuations in the inertial range of the Navier-Stokes case 
(c). The gamma law 10 -3 exp(— \5E\/0A)\5E\~ 3A is represented by the dashed curve. 
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FIG. 5: The PDFs of Elsasser field fluctuations 5z + for the same five different scales as in Fig. [TJ 

the rescaled PDFs according to Eq. ([3]) for MHD case (a) (similar for (b), not shown) while 
Fig. H] displays the rescaled PDFs obtained from the Navier-Stokes simulation (c). The 
corresponding increment distances £ are all lying in the respective inertial range. Evidently 
the PDFs are self-similar and collapse for up to 20a with weak scattering on the master 
PDF, P s , when using the characteristic exponents given above. The dashed lines in both 
figure display the best fitting gamma laws. 

The PDFs of the Elsasser field fluctuations, P[<5,2 + (£)], in system (a) (systems (b) and 
(c) likewise) display a different and well-known behavior as can be seen from Fig. The 
distributions lose their small-scale leptocurtic character as £ increases. Due to the lacking 
correlation of distant turbulent fluctuations the associated distributions become approxi- 



mately Gaussian at large scales. Because of the resulting multifractal scaling of the PDFs 
which is a signature of the intermittent small-scale structure of turbulence it is obvious that 
they can not be collapsed onto a single curve even in the inertial range. However, one can 
infer the nonintermittent characteristic scaling exponent by regarding the function P(0, £) 
(not shown). For example, in system (a) this function exhibits clear inertial-range scaling 
~ £~ a with a = 0.33 ± 1.5 x 10~ 2 in very good agreement with «K4i- 

The occurrence of gamma PDFs is made plausible by a simple reaction-rate ansatz [14j, 



151 ]: Consider the 'intensity' n(e) of turbulent fluctuations with energy e = \SE\ such that 
n(e) is the fraction the total turbulent energy associated with these fluctuations and the 
larger eddies in which they are embedded. The evolution of this function is assumed to obey 
the following linear rate equation: 

de'n(e')/T+(e', e) (4) 

where r + (e',e) is the time characteristic of the creation of fluctuations with energy e as a 
result of turbulent transfer from fluctuations with energy e' while r_(e) is the respective 
characteristic decay time. Normalization of n(e) by J„ de'n(e') yields the corresponding 
PDF -P(e). In a statistically stationary state Eq. (jll) then gives 

P(e) = C i r<k>P(er)-±¥L (5) 

where C\ is a normalization constant. For r_(e)/r + (e', e) ~ (e'/e) 7 this integral equation has 
the solution P(e) = C2e~ 7 exp(— e/A). Thus, the model (jl]) which mimics in combination 
with the abovementioned assumptions a direct spectral transfer process yields the observed 
gamma distributions. Note that the lower bound of the integral in Eq. ([5]) implies that 
energy flows from higher to lower levels where for technical simplicity very large differences 
between e and e' are allowed. A finite upper bound of the integral in Eq. ([5]) does, however, 
not change the result fundamentally. This suggests that the observed gamma distributions 
are an indication of turbulent spectral transfer. 

In summary it has been shown by high-resolution direct numerical simulations of incom- 
pressible turbulent magnetohydrodynamic and Navier-Stokes flows that the monoscaling of 
energy fluctutation PDFs observed in the solar wind is the consequence of lacking Galilei 
invariance of energy increments in combination with self-similar scaling of the underlying 
turbulent fields. The closeness of the PDFs to Levy-type gamma distributions is made 
plausible by a simple model mimicking nonlinear spectral transfer. 
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